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We develop a free energy framework to describe the response of glasses to applied stress. Unlike 
crystals, for which the free energy increases quadratically with strain due to affine displacements, 
for glasses, the nonequilibrium free energy decreases due to complex interplay of non-affine displace¬ 
ments and dissipation. We measure this free energy directly in strained colloidal glasses, and use 
mean-field theory to relate it to affine and nonaffine displacements. Nonaffine displacements grow 
with applied shear due to shear-induced loss of structural connectivity. Our mean-field model allows 
for the first time to disentangle the complex contributions of affine and nonaffine displacements and 
dissipation in the transient deformation of glasses. 


Unlike crystals that deform reversibly under small ap¬ 
plied strain, amorphous materials flow irreversibly even 
at vanishing applied deformation. This generic relaxation 
of glasses is due to complex interplay of affine and non¬ 
affine atomic displacements and dissipation, which are di¬ 
rect consequences of the lack of lattice symmetry. While 
the resulting transient behavior is central to all mechan¬ 
ical properties such as elasticity, time-dependent anelas- 
ticity and flow of glasses, its understanding remains a ma¬ 
jor challenge. An appealing physical picture of relaxation 
is provided by the free energy landscape of glasses. This 
concept has remained, however, rather abstract, lacking 
the direct connection to the atomic displacements. This 
dilemma is fundamentally connected to the lack of lat¬ 
tice symmetry: in crystals, the lattice point symmetry 
guarantees local balance of forces by affine displacements, 
while in glasses, the short-range order precludes force bal¬ 
ance by affine displacements alone, leading to additional 
non-affine motion to restore force balance locally. Dis¬ 
entangling the complex contributions of affine and non¬ 
affine displacements and connecting it to the free energy 
landscape of amorphous materials remains a crucial chal¬ 
lenge. 

Colloidal glasses allow direct observation of single par¬ 
ticle dynamics, and stand out as good model systems 
to study the microscopic degrees of freedom of amor¬ 
phous solids. The particles exhibit dynamic arrest due to 
crowding at volume fractions larger than (fg ^ 0.58, the 
colloidal glass transition [ 1 ]. Recent theoretical [ 2 ] and 
experimental work shows that in these hard-sphere 
glasses, the particle configurations provide a unique route 
to the free energy determined by geometry only. This 
opens up new opportunities to elucidate, experimentally, 
the relation between microscopic degrees of freedom and 
important bulk thermodynamic quantities of deforma¬ 
tion. 

In addition, for the bulk elastic energy, a recent the¬ 
oretical framework generalizes Born’s stability criterion 
from crystals to amorphous solids. The key ingredient 
is the lack of a local center of symmetry in amorphous 
solids, as opposed to crystals: strain-induced forces that 
do balance by symmetry in centrosymmetric crystals, re¬ 


quire additional non-affine displacements in amorphous 
solids to be relaxed. The resulting framework relates 
affine and non-affine displacements to the stability of 
amorphous solids, offering new opportunities to describe 
their transient deformation. 

In this Letter, we provide a generic framework of the 
transient deformation of glasses by connecting micro¬ 
scopic degrees of freedom to the free energy of deforma¬ 
tion. We combine direct measurements of the free energy 
in sheared hard-sphere glasses with the extended Born’s 
stability framework to link affine and nonaffine displace¬ 
ments to the moduli, stress and free energy of the glass. 
By applying this formalism to the yielding of colloidal 
glasses we elucidate how elasticity, vibrational entropy 
and dissipation interrelate in the deformation. We find 
that, contrary to the quadratic free energy increase of lin¬ 
ear elastic solids, for glasses, the non-equilibrium free en¬ 
ergy of deformation decreases with applied strain due to 
elastic energy being lost in nonaffine displacements and 
viscous dissipation. We provide a simple particle-scale 
picture of the transient deformation: the loss of inter¬ 
particle contacts in the extensional sectors of the shear 
plane leads to loss of connectivity and proliferation of 
non-affine displacements, resulting in a marginally stable 
state when the material eventually yields. This frame¬ 
work allows us to disentangle the complex interplay of 
rigidity, dissipation, and disorder in the deformation of 
glasses. 

We use a colloidal glass consisting of sterically stabi¬ 
lized fluorescent polymethylmethacrylate (PMMA) par¬ 
ticles with a diameter of d = 1.5/im, and a polydis- 
perity of 7%, suspended in a density and refractive in¬ 
dex matching mixture of Cycloheptyl Bromide and Cis- 
Decalin. A dense glassy suspension with particle vol¬ 
ume fraction / ^ 59% is prepared by diluting suspen¬ 
sions centrifuged to a sediment, resulting in a relax¬ 
ation time of the colloid of r = 2 x [5]. We 

shear the suspension in between two parallel plates 65/im 
apart at constant rates between 7 = 1.5 x and 

7 = 2.8 X of the order of the inverse relaxation 

time r Individual particles are imaged using confocal 
microscopy and their centers are located in three dimen- 
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FIG. 1 . (Color online) Free energy and volume fraction upon 
deformation of crystal and glass. | (a) [ Sketches of sheared crys¬ 
tal and glasses at ze ro strain (dash boxes) and strain 7 ^ 0.18 
(solid boxes). |(b)| Free energy per particle as a function of 
strain: constructed affine deformation of fee crystal (stars and 
up-curved parabolic fit), and experimental deformation of a 
colloidal glass (filled symbols and down-curved fit according to 
our model). Symbols indicate shear rates of 7 = 1 . 5 xl 0 “^s“^ 
(black squares), 3 x (red circles), 6x (green 

triangles), (blue triangles), 2.8 x (orange 

diamonds). Inset: a sk etch of Voronoi volume (^’^;) and free 
volume (vf) definition. |(c)| Particle volume fraction changes 
Acj) with strains of the colloidal glass at shear rate 
Inset: Theoretical volume fraction changes with strains in a 
simple linear expansion (blue). 


sions with an accuracy of 0.03 /im in the horizontal, and 
0.05 iim in the vertical direction [5]. To measure the 
free energy 012 ], we generate a Voronoi tessellation that 
divides the space into distinct, non-overlapping convex 
polyhedra. We then define the free volume Vfi of parti¬ 
cle i as the volume of a smaller cell generated from the 
Voronoi cell by moving the faces normally inside over a 
distance D = d/2. The total free energy of the hard- 
sphere system can be expressed directly in terms of only 
the free volume as F where 


A = (k?/2i^mkBTY^^ is the thermal wavelength. Note 
that this method is applicable at high densities [2] espe¬ 
cially for glassy system [3] where the particles are densely 
packed. 

The free energy links the microscopic degrees of free¬ 
dom directly to the bulk elasticity: For a uniformly 
sheared crystal, the free energy increases quadratically 
with applied strain, similar to the elastic deformation of a 
spring. This is indeed what we find when we construct an 
affine shear deformation on a computer-gene rated hard- 
sphere face-center-cubic (fee) crystal (see fig. |l(a)[ left), 
and use our method above to compute the free energy 
(see open stars and blue fit in Fig. |l(b)| ). Here, the 
particle volume fraction 0 = 54.5% was preserved us¬ 
ing periodic boundary conditions. The curvature of the 
parabola indicates the elastic modulus, which we deter¬ 
mine to be /i M.2(f/kBT, in good agreement with 
/i ^ 32.4:d^/kBT determined for 0 ^ 54.3% in simula¬ 
tions [6]. In contrast, when we shear a colloidal glass, we 
measure a surprising decrease of the free energy (filled 
symbols and black fit in Fig. |l(b)| ). This monotonic 
decrease indicates that microscopic degrees of freedom 
other than affine dominate the deformation. 

A possible mechanism decreasing the free energy is 
dilation; glasses can dilate under shear, increasing the 
free volume, and hence decreasing the free energy. We 
thus monitored, by particle counting, the volume fraction 
upon deformation; we find, however, that the difference 
A0 to the quiescent state is small and largely constant 
in the strain window where the free energy decreases, 
see Fig. |l(c)| (main panel and inset, black line) [7]. This 
is confirmed by independent check of the corresponding 
pair correlation functions, whose perfect overlap also in¬ 
dicates a closely constant volume fraction, see [8]. In 
contrast, by simple linear expansion for marginal elastic 
solids [ 9 ] we estimate a much stronger dilatancy by as 
muc h as A0 = 40% within the same strain window, see 
Fig. |l(c)| (inset, blue line) and ref. [8]. This linear expres¬ 
sion is defined in terms of a dilatant strain e = \Rp"i^ 
where Rp = -% = (%■ %)i “ f is the Reynolds 

coefficient, E = 2G(1 + v) is the Young modulus, G the 
shear modulus and = | the Possion ratio. The resulting 
strain-dependent volume fraction under constant pres¬ 
sure, (j) = 00(1 + with 00 the volume fraction 

of the undeformed state, decreases considerably, in co n- 
trast to what is observed in the experiment (Fig. 1(c)). 

Hence, the linear expansion grossly overestimates the 
dilation for all but the smallest deformations and we 
conclude that the deformation becomes nonlinear very 
rapidly. Indeed, the distribution of free volumes suggest 
significant redistribution; such redistributions are caused 
by non-affine displacements, since affine displacements 
would leave the overall distribution unaffected. We there¬ 
fore consider nonaffine displacements to redistribute the 
microscopic volumes and to lower the free energy of defor¬ 
mation [TOl [ 11 ] . The main idea is that if the particles are 
not local centers of symmetry (as is the case in all amor¬ 
phous solids), there is an imbalance of forces on every 
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particle when a deformation is applied, unlike in crystals 
with center of inversion symmetry. This additional net 
force acting on every particle leads to additional non- 
affine motions on top of the affine displacements. Be¬ 
cause the nonaffine displacements perform internal work 
against the applied force field, this results in a net nega¬ 
tive contribution to the free energy, contrary to the affine 
displacements that contribute positively. 

To estimate this quantitatively, we first consider the 
affine part of the free energy. Fa = where 

the shear modulus Ga = in the linear regime, 

according to the Born-Huang theory of lattice dynam¬ 
ics [12]. Here, is the number of nearest neighbors and 
K the spring constant associated with a nearest-neighbor 
bond m- The number of bonded neighbors is given by 
the integral of the first peak of g{r), which yields ^ 12 
for the static hard-sphere glass m, but beco mes lower 
under applied shear, as illustrated in Fig. 2(a)[ Particles 
become crowded in the compression sector of the shear 
plane, whereas they become dilated in the extension sec¬ 
tor. Because of strong excluded volume interactions (the 
nearest neighbors cannot come closer to a selected par¬ 
ticle than its excluded volume), the particle increase in 
the compression sector is comparatively small, and does 
not completely balance the particle loss in the extension 
sector, leading to a net loss of particles. 

This is indeed what we observe in the experimental 
pair correlation function reso lved a long the extension and 
compression direction (Fig. |2(b)[ ): The first maximum 
g{r)max decreases in the extension direction, indicating 
particle loss, while in the compression direction, it in¬ 
creases only slightly to saturation. Assuming that the lo¬ 
cal cage dynamics is governed by the Smoluchowski equa¬ 
tion with shear m, we find that the number of nearest 
neighbors (proportional to the peak of g{r)) decreases 
exponentially with strain, 715(7) = 7i^ exp(—Ay) corre¬ 
sponding to a decrease of g{r) in the extension direction 
gd'^) = 3.07exp(—Ay) (see [8]), in excellent agreement 
with the exponential dependence shown in Fig. |2(b)] The 
numerical factor A follows from fit to the experimental 
data for g{r)^ and we obtain A = 0.9. This shear-induced 
cage-breakdown may be even bigger than what we can 
here infer on the basis of the static g{r) [18]. 

The reduced connectivity leads to growing non-affine 
contributions to the shear modulus. In the spirit of 
Alexander m, the nonlinear (non-affine) contribution 
to the shear modulus can be written as resulting from a 
Taylor expansion in the free energy up to third order in 
7, Gna = (6 + C7), where C is a phenomenological 

constant from the non-affine free energy expansion. 

We can now evaluate the free energy of deforma¬ 
tion quantitatively. The total nonequilibrium free en¬ 
ergy of deformation is F(y) = Fq + Fei{j) - Wdissil), 
where Fq is the strain-independent, ’’ground-state” en¬ 
ergy of the metastable minimum (inherent structure) of 
the glass, Wdiss{l) the energy dissipated due to the mi¬ 
croscopic friction and irreversible rearrangements, and 
Fei = FAdi) — Fa[a{f) is the elastic (reversible) en- 



FIG. 2. (Color online) Shear-induced loss of structural con¬ 
nectivity. (a) Schematic of nearest neighbor configuration 
with applied shear. Upon application of shear, the number of 
particles moving out of the cage in the extension direction is 
larger than that moving in along the compression direction, 
leading to net loss of connectivity, (b) Measured maximum 
of g{r) in the compression (black squares) and extension (red 
circles) axis (log scale) as a function of strain (linear scale). 
Black fitting line has slope —0.9. 



FIG. 3. (Color online) Comparison of model prediction (blue 
line) and experimental data (black squares) of stress (cr) as 
a function of strain (log-log scale). To measure the stress 
reliably at these small shear rates, we used smaller particles 
with diameter di = 50nm, resulting in shear stresses larger 
by a factor [d/di ~ 3 X 10^. Data was measured at same 
volume fraction cj) = 0.59 and similar Peclet number. Black 
line has slope 1 to guide the eye. Inset: Theoretical prediction 
of the average number of nearest neighbors (log scale) as a 
function of strain (linear scale). 


ergy, which can be simply expressed by the shear mod¬ 
ulus as Fei = with G = Ga — Gna^ containing 

both affine and non-affine contributions. The dissipated 
energy associated with a continuous ramp of strain at 
strain rate 7 is Wdiss = Jq cr{s)jds^ where the stress 

cr{t) =7 Jq G{s)ds [20] with G the (strain dependent) re¬ 
laxation modulus. Using a standard viscoelastic model, 
the relaxation modulus has the general form G{t) = 
G + Gr exp[—t/r]^, where r = ^/Gr is the global relax¬ 
ation time, 77 the viscosity m and Gr = Gq — G, with 
Go the instantaneous (infinite-frequency) shear modulus. 
For the standard linear viscoelastic solid [20], one has 
/3 = 1, while for many glassy materials the relaxation 
is stretched exponential with /3 < 1. Thus, inserting 
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FIG. 4. (Color online) Minkowski measure of cage-distortion 
underlying the loss of connectivity. Eigenvalue ratio /3i° of the 
Minkowski tensor indicates elliptical cage distortion, (a) 
Average ratio (/3i°) as a function of strain for 7 = 
Decreasing values and plateau indicate increasing shear dis¬ 
tortion and steady state, (b) Steady state values of (/3i°) 
(circles) and (p/q) (squares) as a function of shear rates. De¬ 
creasing values indicate increasing elliptical distortion. Inset: 
Reconstructed effective elliptical cage. Color indicates ap¬ 
plied shear rates 0.6 (black), 1.5 (red), 3 (green), 6 (blue), 10 
(cyan) and 28 (purple) xl 0 “^s“^. 


the stretched-exponential expression for G{t) in the inte¬ 
gral, and using t = 7/7, one obtains Wdiss as a function 
of 7, 7 and Gr. Focusing on the limit of very low 7, 
the leading term in a Taylor expansion around 7 = 0 is 
Wdiss ~ independent of /?. 

Having specified the various contributions, the final 
expression for the nonequilibrium free energy of the col¬ 
loidal glass subject to a strain ramp is 

F(7) = Fo + -^ — [n‘^ exp(-^7) - (6 + 

DTI (J I 

( 1 ) 

For a crystal with center-of-inversion symmetry, the non- 
affine part F/va vanishes and, because typically Ga > 
Gr^ the free energy becomes a monotonic increasing func- 
tion of 7, just like the crystal free energy in Fig. |l(b)| 
(open stars and blue line). In amorphous solids, however, 
non-affine contributions are significant, and ultimately 
lead to the decreasing behavior of the free energy (jointly 
with the effect of dissipative terms, also negative). To 
test our model quantitat ively, we simultaneously fit the 
free energy data in Fig. |l(b)| using Eq.(l) (black line) 
and stress-strain curves of a colloidal glass measured in¬ 
dependently with a rheometer, see Fig. We compute 
the stress as cr = o^i^Gdiss^ where we take = dFeijd^^ 
and Gdiss = 777(1 — 6“^/'^^). We obtain excellent simulta¬ 
neous fits up to large strains, lending credence to the re¬ 


lation between nonequilibrium free energy and non-affine 
and dissipative contributions. The predicted number of 
nearest neighbors (Fig. inset) demonstrates the loss of 
connectivity as a function of strain. 

We finally highlight the structural distortion of the 
strained glass using the Minkowski formalism |22]. Un¬ 
like g{r) that represents a statistically averaged neigh¬ 
bor distribution, the Minkowski metrics measure directly 
the local cage anisotropy from higher moments of the 
Voronoi volume and surface distributions. The tensor 
r^dA, where r is the vector from the central 
particle to its Voronoi neighbors and S is the curvature 
of the infinitesimal area element dA, quantifies elliptical 
distortions: the average ratio of the smallest and 

largest eigenvalues of is directly related to the ratio 
of the two semi-axes p and q of an effective elliptical cage 
according to {Pi^) = [p/q^] a value of {Pi^) = 1 indi¬ 
cates an effective isotropic cage, while {Pi^) < 1 indicates 
elliptical distortion. Indeed, the average of reveals in¬ 
creasing elliptical distortion of the nearest neighbor cage, 
and saturation at high strains indicating steady state (see 
fig- El (a)). The steady-state distortion as a function of 
strain rates, shown in Fig. (b), reveals cross-over into 
shear-rate independent behavior at 7 ^ . The recon¬ 

structed effective elliptical cages (inset) demonstrate the 
increasing nearest-neighbor distortion. The plateau at 
7 ^ is in agreement with our assumption that the 
cage distortion and number of nearest neighbors become 
independent of strain rate in the low strain-rate regime 
investigated here. 

In summary, our novel framework connecting micro¬ 
scopic degrees of freedom to the nonequilibrium free en¬ 
ergy of glasses allows quantitative description of the tran¬ 
sient deformation of glasses far into the nonlinear regime. 
The sum of the three free energy contributions, affine, 
non-affine and dissipation, explains the surprising de¬ 
crease of the nonequilibrium free energy with strain, un¬ 
like in crystals. While our hard-sphere colloidal glass 
allows direct measurement of the free energy and under¬ 
lying microscopic distortions, the proposed mechanism 
should apply to molecular glasses as well. For those 
molecular glasses, the strong hard-core repulsion is re¬ 
placed by the steep repulsive potential (oc r^^) due to 
Pauli’s principle, leading to similar shear-induced loss of 
connectivity and growth of non-affine dynamics. Hence, 
interestingly, it is precisely the quantum mechanical prin¬ 
ciple guaranteeing atomic stability that destabilizes the 
amorphous material on a larger scale. 
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